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To see whether heavy baryon and anti-baryon can form a 
bound state, the heavy baryonium, we study the two-pion ex- 
change interaction potential between them within the heavy 
baryon chiral perturbation theory. The obtained potential is 
applied to calculate the heavy baryonium masses by solving 
the Schrodinger equation. We find it is true that the heavy 
baryonium may exist in a reasonable choice of input param- 
eters. The uncertainties remaining in the potential and their 
influences on the heavy baryonium mass spectrum are dis- 
cussed. 

1 Introduction 

Quark model has achieved great success in describing the experimentally observed 
hadronic structures to a large extent. And the quark potential in between quark and 
anti-quark deduced from Chromodynamics (QCD) can explain the meson spectrum quite 
well. Many of predicted states by potential model were discovered in experiment and 
the theoretical predictions are in good agreement with experimental data, especially in 
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charmonium and bottomonium sectors [H [21 E], where the masses of charm and bot- 
tom quarks are heavy enough to be treated non-relativistically. However, things became 
confused after the discovery of X(3872) in 2003 at Belle [1], which was later confirmed 
by BaBar [5]. In recent years, a series of unusual states in charmonium sector, such 
as F(4260), y(4360), F(4660), and Z±(4430), were observed in experiment [6]. Due to 
their extraordinary decay nature, it is hard to embed them into the conventional char- 
monium spectrum, which leads people to treat them as exotic rather than quark-quark 
bound states. The typical scenarios in explaining these newly found states include treat- 
ing y(4260) as a hybrid charmonium [7], a XcP^ molecular state [8], a conventional \E'(4S') 
[9], an uJXci molecular state [10], a AcAc baryonium state [11], a DiD or DqD* hadronic 
molecule [12], and a P-wave tetraquark [cs][cs] state [I3]; F(4360) is interpreted as the 
candidate of the charmonium hybrid or an excited D-wave charmonium state, the 3^Di 
[Hj and an excited state of baryonium [16]; y(4660) is suggested to be the excited S-wave 
charmonium states, the 5^5*1 [14J and 6^Si \X5\, a baryonium state [16l[T7], a /o(980)\E'' 
bound state [THJ [19], a d^Si-i^Di mixing state [20], and also a tetraquark state [211 [22] . 
There have been recently many research works on "exotic" heavy quarkonium study in 
experiment and theory. To know more of recent progress in this respect and to have a 
more complete list of references one can see e.g. recent reviews [23l |2l] and references 
therein. 

In the baryonium picture, the tri-quark clusters are baryon-like, but not necessar- 
ily colorless. In the pioneer works of heavy baryonium for the interpretation of newly 
observed "exotic" structures [TTl [TB] , there were only phenomenological and kinematic 
analysis, but without dynamics. In this work we attempt to study the heavy baryonium 
interaction potential arising from two-pion exchanges in the framework of Heavy Baryon 
Chiral Perturbation Theory (HBCPT) [23]. The paper is organized as follows. In Section 
2, we present the formalism for the heavy baryon-baryon interaction study; in Section 3 
we perform the numerical study for the mass spectrum of the possible baryonium with 
the obtained potential in preceding section; the Section 4 is devoted to the summary and 
conclusions. For the sake of reader's convenience some of the used formulae are given in 
the Appendix. 
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2 Formalism 



To obtain the heavy baryonium mass spectrum, we first start from extracting the 
baryon-baryon interaction potential in the same procedure as for quark-quark interaction 



2.1 Heavy Baryonium 

In the heavy baryonium picture [16], and S° are taken as basis vectors in two- 
dimensional space. The baryonia are loosely bound states of heavy baryon and anti- 
baryon, namely 

Bt ^ |A+S°> 
Triplet : fi? ^ -^(|A+ A+ > - IS^SO >) (1) 
^ |A+S°> 

and 

Singlet : 5° = ^i\K K > + l^c^c >) • (2) 

Here, approximately the transformation in this two-dimensional "C-spin" space is invari- 
ant, which is in analog to the invar iance of isospin transformation in proton and neutron 
system. 



2.2 Effective Chiral Lagrangian 

Heavy baryon contains both light and heavy quarks, of which the light component ex- 
hibits the chiral property and the heavy component exhibits heavy symmetry. Therefore, 
it is plausible to tackle the problem of heavy baryon interaction through the heavy chiral 
perturbation theory. Following we briefly review the gists of the HBCPT for later use. 

In usual chiral perturbation theory, the nonlinear chiral symmetry is realized by mak- 
ing use of the unitary matrix 



where M is a 3 x 3 matrix composed of eight Goldstone-boson fields, i.e., 



M 



V 



72" "^76^ 

K- 



TT 

. J_7rO I J_„ 



(4) 



Here, /^r is the pzon decay constant. 



After the chiral symmetry spontaneously broken, the Goldstone boson interaction with 
hadron is introduced through a new matrix [26l [27] 

(5) 



1 iM_ 

^ = S 2 = e 



From ^ one can construct a vector field and an axial vector field with simple chiral 



transformation properties, i.e. 



(6) 
(7) 



For our aim, we work only on the leading order vector and axial vector fields in the 
expansion of ^ in terms of they are 

1 



= —Md^M , 

J TV 



A,, 



In 



(8) 
(9) 



For heavy baryon, each of the two light quarks is in a triplet of fiavor SU(3), and hence 
the baryons can be grouped in two different SU(3) multiplets, the sixtet and antitriplet. 
The symmetric sixtet and antisymmetric triplet can be constructed out in 3 x 3 matrices 
they are 



/ V++ J_y+ J_"'+ \ 
J_y+ yO 1 ■^'0 

\ ^2 c c c 



(10) 



and 






Ae 


' — 'c 


-Ac 





' — 'c 


' — 'c 


' — 'c 






respectively. 

Introducing six coupling constant Qi, i = 1,6, the general chiral-invariant Lagrangian 
then reads 



Cg = ^tr[B-,{tI^-M-^)B-,]+tr[BeitI^-Me)Be] 

+ tr[BZ[-g^,{iI) - M*) + i{^^D, + ^,D^) - 7^(^i? + M*)7.]5r] 

+ qMB^Ip^I^^^Bq) + g2tr{B^lf,l^A^'B-^) + h.c. 

+ gMBl^A^B^) + h.c. + gMBl^A^B-^) + h.c. 

+ (75tr(S6"*7M75^^56J+^?6tr(537^75^'^53) • (12) 

Here, B^^ is a Rarita-Schwinger vector-spinor field for spin-| particle; M3, Mg, Mg rep- 
resent for heavy baryon mass matrices of corresponding fields; With the help of vector 
current defined in Eq. (jH]), we may construct the covariant derivative D^, which acts 
on baryon field, as 

D^B^ = d^B, + V^B^ + B^V^ , (13) 

D,B-^ = d,B-, + V,B-^ + B-^V^ , (14) 

where stands for the transpose of V^. Thus, the couplings of vector current to heavy 
baryons relevant to our task take the following form 

Cs, = \tr{B-,tYV,B-,) 

^ Ae27^(7r°9^7r° + TX-d^Ti+ + 7r+9^7r-)Ae , (15) 



and 



le, = UriB-.B-.iYV;^ 
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= ^K^K^t^'^i^n^dy + Ti-d^ + TT+S^TT-) . (16) 

According to the heavy quark symmetry, there are four constraint relations among those 
six coupling constants of the Lagrangian of Eq. (fT2|) . i.e., 

^6 = , ^3 = , 95 = -^91 , 9i = -V3g2 , (17) 
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Figure 1: Schematic Diagrams which contribute to the baryonium potential 



which means the number of independent couphngs are then reduced to two. In this work, 
we employ gi and g2 for the numerical evaluation as did in Ref. |25j . 

Here, to get the dominant interaction potential we restrict our effort only on the pion 
exchange processes as usual. Notice that the couplings of pion to spin-| and -| baryons, 
and pion to two spin-| baryons take a similar form, in the following we merely present 
the spin-| and -| harjon-pion coupling for illustration, i.e.. 



£i = °*^9.7r°S° + h.c. 



£2 = +*'^9^7r+S° + h.c. , (19) 

£3 = ++*^9^7r+A+ + h.c. , (20) 
A = -^S °*'^9^7r-A+ + h.c. , (21) 

£5 = -^S +*'^9^7r°A+ + h.c. . (22) 

To get the pion and two spin-| baryon couplings one only needs to replace the S*'* by Sc, 
by by (72? and insert 7'^75 in between the two baryon fields in Eqs.f lT8|) - fl22|) . 



2.3 Baryonium Potential from Two-pion Exchange 

To obtain heavy baryon-baryon interaction potential in configuration space, we start 
from writing down the two-body scattering amplitude in the center-of-mass frame(CMS), 
i.e. taking = — pj, and p'„ = — p^. In CMS the total and relative four momenta are 
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defined as 



P = {Pa + Pb) = (P'a + P'b) = (E, 0) , (23) 

P = liPa - Pb) = (0, p) , (24) 

P' = lip'a - p'b) = (0, P) . (25) 

To perform tlie calculation, it is convenient to introduce some new variables as functions 
of p and p', i.e., 

W(p) =EM + Eb{p), (26) 
W(p') = Ea{p') + E,{p') , (27) 
FEip,Po) =^E+po-Eip)+t5 , (28) 

where 6 is an infinitesimal quantity introduced in the so-called i6 prescription. Following 
the same procedure as in Refs. [281 129], it is straightforward to write down the baryon- 
baryon scattering kernels, shown as box and crossed diagrams in Figure [H 

Kto.= - j^{E-W{p')){E-W{p)) J dp'odpodk2odhod'k2d'k 



X 



:S^{p-p -h- k2) 



1 

1 1 



(2vr)4 ^ ^>kl-m^ + i5FE{v',p',)FE{-v',-p'o) 



Fe{p - k,j9o - A;io)F£;(-p + k, -p^ + fcio) Fe(jp,Pq)Fe{v, -po)) 
^ r2 — W^' (29) 

Kcross= - -^{E-W{p')){E-W{p)) J dp'odpodk2odkiod^k2d^ki 



X 



X 



rS'^ip-p - ki- k2 



(2vr)4 ^ ''kl-m^ + z6FE{p',p'o)FE{-p',-p'o) 

Fe(p - k,po - A;io)Fi5(-p' - k, -p'^ - kio) Fij(p,po)^£;(-p, -Po) 
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Here, m corresponds to the pion mass and r^j- are heavy haijon-pion interaction vertices 
that can be read out from the Lagrangian in Eqs. (ll8p - (l22l) . In case of spin-| intermediate, 

TjTiTiTj = (^Y^ u{-p)k!^u^{p - ki)ui,{p - ki)k'[u{p) 

X vip)i-k'^)vai-p + ki)vpi-p + k,)klv{-p) , (31) 
and in case of spin-| intermediate 

TjTiTiTj = (^j-^ u{-p)-ff,-f5k!^u{p - ki)u{p - fci ) 7^.75 /cim(p) 

X v{p)-fo,-f^{-k'l)v{-p + ki)v{-p + /ci)7/375^2^(-p) • (32) 

Integrating over p'q, pq, kiQ, and /C20 in Eq.f E^ one obtains the interaction kernel of 
box diagram at order C(-p^), 

1 f d^kid^k2 TjTi 



box 



(27r)3 J AE^.E^, Ep_kx + Ep-W + E^, 

E^ + Ep_k, -W + E^,Ep + Ep,-W + E^^ + ' ^^^^ 
where Mh represents one of the heavy baryon mass, M^+, M^o or Ms*; -Ep-ki = 



p — ki)2 + M|,, is the intermediate state energy; i?ki = V^i + -^k2 = 

are two pions' energies; and = 2E{-p). With the same procedure, we can get the in- 
teraction kernel of crossed diagram, i.e. 

1 f d^kid^\i2 TjTi 



Across ^27r)3 7 4Ek,£;k. ^p-k, + ^p-Vr + ^ki 



r r 1 

(34) 



E' + Ep,+ki - W + E^,E^ + E^.-W + Ek, + 



Next, since what we are interested in is the heavy baryons, we can further implement 
the non-relativist ic reduction on spinors with the help of vertices given in Eqs. (fT8|) - (l22|) . 
In the end, the non-relativistic reduction for A^S+*7r'^ and A^S^tTr" couplings gives 

^ (f^) ^(^2)?^m(Pi)(P2 - Pi)'' = -i S'^ ■ q , (35) 

(^j-^ m(p2)7m75m(Pi)(P2 - ViY = i (^j-^ o-i ■ q , (36) 



and 

i 



respectively. Here, q = P2 — Pi and is the spin-i to spin-| transition operator. 

In the process of deriving — A+ potential, the S+ and S+* are taken into account 
as intermediate states. Using Eqs. (!35|) - (l36|) and the explicit forms of spinors given in the 
appendix, we can readily obtain the reduction forms for the S+ intermediate 

u{-p)l^il5Ku{p - ki)u{p- ki)-fi,-f5k'[u{p) X 
viphal5i-ki)v{-p + ki)v{-p + /ci) 7/375 
= (ki-k2)' + (cri.ki X k2)(cr2-ki X ka) , (37) 

the E^* intermediate in the box diagram 

u{-p)k(^u^{p - ki)u^{p - ki)k\u{p) x 
v{p){-k'^)va{-p + ki)vp{-p + ki)k^v{-p) 
= ^(ki-k2)2-^(o-i-ki X k2)(o-2-ki X ka) , (38) 

and the crossed diagram 

u{-p)k^Uf,{p - ki)u^{p - ki)klu{p) x 
v{p){-k'^)va{-p + ki)vp{-p + ki)k^v{-p) 

= ^(ki-k2)2 + i(<Ti-ki X k2)(0-2-ki X k2) , (39) 

respectively. Thus, the spinor reduction finally leads to an operator Oi(ki, k2), of which 
the variables ki and k2 can be replaced by gradient operators Vi and V2 in configuration 
space and acting on ri and r2, respectively. This operator is expressed as 

Oi(ki, k2) = CiOi(ki, k2)+C202(ki, k2) 

= ci(ki ■ k2)^ + C2(cri • ki X k2)(cr2 • ki X ka) . (40) 

Here, the decomposition coefficients Ci and C2 are given in Table 1. The first part of 
Eq. fHOj) may generate the central potential and the second part may generate the spin- 
spin coupling and the tensor potentials, which are explicitly shown in the Appendix. 
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Table 1: The values of coefficients ci and C2 in the decomposition of operator 0(ki, in 
Eq. (j40p . The left one is for the spin-^ intermediate state case and the right one is for the spin-| 
case. 



spin- 1/2 ci C2 


spin-3/2 ci C2 


box 1 1 
cross 1 1 


box 4/9 -1/9 
cross 4/9 1/9 



To get the leading order central potential, e.g. for Ac-Ac system, we first expand the 
energy in powers of but keep only the leading term, like 

I « I ^ ^ (41) 

where Ai = M^* — Ma^ represents the mass splitting. By virtue of the factorization in 
integrals given in the Appendix, we can then make a double Fourier transformation, i.e.. 



Vcf (ri, r2) 



9t 




(42) 



/4; J J (27r)6 2E^,E^,{E^, + Ai){Eu, + Ai){E^, + E^, 

where the superscript B denotes the box diagram and the subscript C means central po- 
tential. Similarly, one can get the central potential from the crossed diagram contribution 

J^J J J ^2^ 



V^in, r2) 



Oi(ki,k2)e^'^i'-^e^'^^''^/(k^)/(k^)D 



(43) 



where the superscript C denote crossed diagram and the subscript C means central po- 
tential, and 



D 



1 



+ 



4-Eki -£^k2 

1 



1 



+ 



1 



(Ek, + Ai)2 (Ek, + Ai)V + ^k, 
1 2 



1 



+ 



+ 



(44) 



,(^ki + Ai)2 (Ek, + Ai)2 (Ek, + Ai)(Ek, + Ai); Ek, + Sk2 + 2AiJ 

In order to regulate the potentials we have introduced form factors at each baryon-pion 
vertex. The resulting /(k^) form factors appearing in Eqs. ( 142|) and (143|) will be given in 
Section 3. 

Taking a similar approach as given in above one can readily get the central potential 
in other interaction channels and also the tensor potential. Notice that although there 
exists the one-pion exchange contribution in Sc-Sc system, due to the 7^75 nature in 
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Figure 2: The triangle and two-pion loop diagrams. 

interaction vertex, it only contributes to cri ■ (T2 term, which is out of our concern in this 
work. Here we just focus on the central potential. 

Besides box and crossed diagrams, there are also contributions from triangle and 
two-pion loop diagrams as shown in Fig. [2l As in the box and crossed diagrams, after 
integrating over energy component, we get the pion-pair contribution, as shown in the 
left diagram of Figure |2l as [33] 



where the 02(^.1,^2) = (ki • k2) from spinor reduction can be replaced in configuration 
space by the gradient operator (Vi ■ V2). Similarly, the two-pion loop contribution, as 
shown in the right diagram of Figure [2] reads 



V2.-ioopirur2) = -[^4 J J ^^^e-^''^e-^^V(kf)/(k^)A . (46) 

Here, A = —jE^ — e^, +Ek ' Expressing Eps. and fll^ in the integral 

representation of -Eki, and making the Fourier transformation, one can then obtain the 
corresponding potentials. 



3 Numerical Analysis 

With the central potentials obtained in preceding section, one can calculate the heavy 
baryonium spectrum by solving the Schrodinger equation. In our numerical evaluation, 
the Matlab based package Matslise [31] is employed. The following inputs from Particle 
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Data Book [321 are used in the numerical calculation: 



M, 



^+ - 2.286GeV , Mso = 2.454GeV , Ms* = 2.518GeV , U = 0.132GeV , m = 0.135GeV 



(47) 



and both spin-^ and -| fermion intermediates are taken into account. 

It is obvious that the main uncertainties in the evaluation of heavy baryonium remain 
in the couplings of Eq. (fTTj) . The magnitudes of the two independent couplings gi and g2 
were phenomenologically analyzed in Ref. [25], and two choices for them were suggested, 
i.e., 

'' = \ = 

and 

(49) 



(4J 



9i 



ixO.75, ^?2 = -J^x0.75 



which implies the g^^ lies in the scope of 1 to 1.4, similar as estimated by Ref. [30] in the 
chiral limit. 

3.1 Gaussian form factor case 



The central potential from two-pion exchange box which can be regularized by widely 
used Gaussian form factor /(k^) = e^'^^^^^ reads 

i dX 



fA 

2Ai 

TT 



n Jq Af + X 
2 Oi(ki,k2 



2 + A2 
°° dX 



Oi(ki,k2)F(A,ri)F(A,r2) 
°° dX_ 

a^Ta^ 



A? + A2^^^'^^W_ aTXT^^(^'^^) 



(50) 



Details of the derivation of Eq. (!50|) from Eq. (H2|) can be found in the Appendix. There, 
the function F{X,r) is defined by Eq. fl72|) . And, similarly the central potential from 
two-pion exchange crossed diagram gives 

1 rdX{Al-X^) 



-Oi(ki,k2)F(A,ri)F(A,r2) 



(51) 
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Here, the ellipsis represents the high singular terms in r2 — )■ ri = r limit, which behave as 
higher order corrections to the potential and will not be taken into account in this work, 
but will be discussed elsewhere. The central potential of Eq. ( 150|) is obtained in the case 
of spin- 1 intermediate state, and the explicit forms of Vcci from box diagram are 

V(fc = p— e-— , 52 

V- , = ^ e-^ , (53) 

32y27r3/^Afr3/2 

16v^7r3/4A?r5/2^ 1287r5/2/4A2r5/2^ • ^^^^ 

With Gaussian form factors it is seen from Eq. fl72l) in the Appendix that for a given A 
the function F(A,r) is suppressed for large A values, that is the dominant contribution 
to potential comes from the small A region. So, in obtaining the analytic expressions 
of above potentials and hereafter, we expand the corresponding functions, as defined 
in the Appendix, in A and keep only the leading term. In this approach, the crossed 
diagram contributes to the potential the same as the box diagram at the leading order in 
A expansion, and hence is not presented here. 

Similarly, we obtain the potentials from triangle and two-pion loop diagrams, i.e.. 



VcGb — ^ ~^ ^ — —e 4 



32y27r7/2/4A^r2^ ' ~ IGV^n^f^Air^/^^ 



1287r5/2/4A,r5/2^ ' ^^""^ 



and 



m^/^A^ 1.2^2 



32y27r2/4r3/2 ^ ' 



To get the central potential for the case of spin-| intermediate state, one needs only 
to make the following replacement 

(?4 ^ (72 , Ai ^ a; = Ms, - Ma, (58) 
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r(fm) r(fm) 

Figure 3: The Ac-Ac central potential behavior in case of Gaussian form factor versus 
different parameter choices. 



in Eq.dSO]). 

Note that in above asymptotic expressions we keep only those terms up to order 
and more singular terms are not taken into accounted in this work. The dependence of 
potential with various parameters are shown in Figure [31 The results indicate that the 
potential approaches to zero quickly in long range in every case, while in short range the 
potential diverges very much with different parameters, as expected. As a result, the 
binding energy heavily depends on input parameters, the coupling constants and cutoff. 
One can read from the figure that in the small coupling situation, the potential becomes 
too narrow and shallow to bind two heavy baryons. Table [2] presents the binding energies 
of Ac-Ac and Sc-Sc systems with different inputs. Schematically, the radial wave functions 
for the ground state of Ac-Ac system with Gaussian and monopole form factors are shown 
in Figure m respectively, while the wave functions for Sc-Sc system exhibit similar curves. 



3.2 Monopole form factor case 

In order to regulate the singularities at the origin in configuration space, usually people 
employ three types form factors in the literature, i.e. the Gaussian, the monopole, and the 
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Table 2: Binding energies with different inputs witli Gaussian form factor. The left table is for 
the Ac-Ac system, and the right one for Sc-^c system. 



1^2! 


A(GeV) 


Binding 


Baryonium 


91 


A(GeV) 


Binding 


Baryonium 






energy 


mass 






energy 


mass 


<0.9 


<0.6 


No 




< 1.0 


< 0.8 


No 




0.9 


0.6 


-22 MeV 


4.550 GeV 


1.0 


0.8 


-11 MeV 


4.895 GeV 


0.95 


0.6 


-77 MeV 


4.495 GeV 


1.05 


0.8 


-61 MeV 


4.845 GeV 


1.0 


0.6 


-168 MeV 


4.404 GeV 


1.1 


0.8 


-145 MeV 


4.761 GeV 


0.95 


0.7 


-196 MeV 


4.376 GeV 


1.05 


0.85 


-141 MeV 


4.765 GeV 


0.95 


0.8 


-227 MeV 


4.345 GeV 


1.05 


0.9 


-266 MeV 


4.640 GeV 


0.95 


0.9 


-588 MeV 


3.984 GeV 


1.05 


0.95 


-438 MeV 


4.468 GeV 




r (fm) r (fm) 

Figure 4: Radial wave function of Ac-Ac ground state. The left figure is for case of 
Gaussian form factor under the condition of \g2\ = 0.95 and A = 0.8, and the right one is 
for the case of monopole form factor with \g2\ = 0.9 and A = 0.95. 



dipole form factors [34]. For comparison we also calculate the potential with monopole 
form factor using the same factorization technique, and the basic Fourier transformation 
for monopole form factor is presented in Appendix for the sake of convenience. Here, in 
obtaining the analytic expressions for potentials we also take the measure of expanding 
the corresponding functions in parameter A and keeping only the leading term. Then, 
what obtained from the box-diagram contribution reads 

T/S _ _ 94 / 'JJ^ -2mr , m ^Ar 

^CM[r)- 8^5/2^4^2^5/2^ 4 ^ ^ 4 

8y27r5/2/^V^^rrAA2r5/2 
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-gj = 0.7A = 0.9GeV 
g2 = 0.7A = 0.95GeV 
g =0.7A = 1.0GeV 



0.8 1.0 

r(fm) 



1.2 



1.4 



> 




Figure 5: The Ac-Ac central potential behavior in case of monopole form factor versus 
different choices of inputs. 



Contributions from triangle and two-pion loop diagrams are 
and 



(60) 



VUr) = -^^lf^IpI!^e~^--^> + iA;zZ^e-- (61) 

respectively, where superscript B, T, and L stand for box, triangle and 2tt loop. Note 
that since there is no heavy baryon intermediate state in the 2tt loop process, as shown 
in the right graph of Figure 2, the potential range of it appears different. 

We find that the structure of potential with monopole form factor is much simpler 
than the Gaussian case. The dependence of potential with various parameters are shown 
in FiglHl From the figure one can see that in small coupling case the potential change less, 
which means the potential tends to be insensitive to the small coupling, and hence the 
binding energy. Solving the Schrodinger equation we then obtain eigenvalues for different 
input parameters, given in Table [31 From the table, we notice that the binding energy 
is sensitive to and changes greatly with the variation of gi, \g2\ and the cutoff A, the 
same as the case with Gaussian form factor. Intuitively, the realistic baryonium can only 
accommodate small ones of those parameters. 
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Table 3: Binding energies with different inputs with monopole form factor. The left table is for 
the Ac-Ac system, and the right one for Tic-Tic system. 



1^2! 


A(GeV) 


Binding 


Baryonium 


91 


A(GeV) 


Binding 


Baryonium 






energy 


mass 






energy 


mass 


<0.7 


<0.9 


No 




< 0.9 


< 0.9 


No 




0.8 


0.95 


-117 MeV 


4.455 GeV 


0.95 


0.95 


-438 MeV 


4.468 GeV 


0.85 


0.95 


-420 MeV 


4.152 GeV 


1.0 


0.95 


-830 MeV 


4.076 GeV 


0.9 


0.95 


-521 MeV 


4.051 GeV 


1.05 


0.95 


-1003 MeV 


3.903 GeV 


0.7 


0.9 


-5 MeV 


4.567 GeV 


0.9 


0.9 


-40 MeV 


4.866 GeV 


0.7 


0.95 


-67 MeV 


4.505 GeV 


0.9 


0.95 


-153 MeV 


4.753 GeV 


0.7 


1.0 


-252 MeV 


4.320 GeV 


0.9 


1.0 


-345 MeV 


4.561 GeV 



3.3 Ground state of A^j-A^ baryonium 



Table 4: Binding energies with the change of parameters for A^-Afe system. The left table is for 
the Gaussian form factor, and the right one for the monopole form factor. Here gi, corresponds 
to g2 in charmed baryonium sector 



\9b\ 


A(GeV) 


binding 


Baryonium 


\9b\ 


A(MeV) 


Binding 


Baryonium 






energy 


mass 






energy 


mass 


<0.7 


<0.7 


No 


No 


< 1.0 


< 0.8 


No 


No 


0.7 


0.75 


-4 MeV 


11.236 GeV 


1.0 


0.8 


-11 MeV 


11.229 GeV 


0.8 


0.75 


-76 MeV 


11.164 GeV 


1.05 


0.8 


-56 Mev 


11.184 GeV 


0.9 


0.75 


-294 MeV 


10.946 GeV 


1.1 


0.8 


-143 MeV 


11.097 GeV 


0.8 


0.8 


-164 MeV 


11.706 GeV 


1.05 


0.8 


-103 Mev 


11.137 GeV 


0.8 


0.9 


-396 MeV 


10.844 GeV 


1.05 


0.9 


-164 MeV 


11.076 GeV 


0.8 


1.0 


-622 MeV 


10.618 GeV 


1.05 


1.0 


-321 MeV 


10.919 GeV 



We also estimate the ground state of Af,-Ah baryonium system with Gaussian and 
monopole form factors. The result are shown in Table HI where Qf, corresponds to g2 in 
charmed baryonium sector. Note that since the dominant decay mode of is to Afovr, 
by which we may constrain the SfeA^vr coupling from the experiment result, and this may 
shed lights on the further investigation on the nature of possible baryonium. 
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4 Summary and Conclusions 



In the framework of heavy baryon chiral perturbation theory we have studied the heavy 
baryon-baryon interaction, and obtained the interaction potential, the central potential, 
in the case of two-pion exchange. The Gaussian and monopole types form factors are 
employed to regularized the loop integrals in the calculation. As a leading order analysis, 
the tensor potential and higher order contributions in expansion are neglected. As 
expected, we found that the potential is sensitive to the baryon-pion couplings and the 
energy cutoff A used in the form factor. 

We apply the obtained potential to the Schrodinger equation in attempting to see 
whether the attraction of two-pion-exchange potential is large enough to constrain two 
heavy baryons into a baryonium. We find it true for a reasonable choice of cutoff A and 
baryon-pion couplings, which is quite different from the conclusion of a recent work in the 
study of DD potential through two-pion exchange |3S] ■ Since usually the cutoff A is taken 
to be less than the nucleon mass, i.e. about 1 GeV in the literature, in our calculation 
we adopt a similar value employed in the nucleon-nucleon case. In Ref. [53] authors took 
a fixed coupling g = 0.59 and obtained the binding with a large cutoff. While in our 
calculation for the baryonium system with Gaussian form factor, there will be no binding 
in case gi < 1.0 and A < 0.8. The increase of coupling constant will lead to an even 
smaller A for a given binding energy. 

Based on our calculation results it is interesting to note that in case there exists binding 
in Sc-Sc system, with both Gaussian and monopole factors, the coupling gi will be much 
bigger than what conjectured in Ref. [25]. However, for A^-Ac system, to form a bound 
state the baryon- Goldstone coupling g2 could be similar in magnitude as what estimated 
in the literature. 

Notice that the potential depends not only on coupling constants and cutoff A, it 
also depends on the types of form factors employed. Our calculation indicates that the 
Gaussian form factor and Monopole form factor are similar in regulating the singularities 
at origin, and lead to similar results, with only subtle difference, for both Ac and Af, 
systems. Numerical result tells that the heavy baryon-baryon potentials are more sensitive 
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to the coupling constants in the case of Monopole form factor, but more sensitive to the 
cutoff A in the case of Gaussian form factor. From our calculation it is tempting to 
conjecture that the recently observed states y(4260) and F(4360), but not F(4660) [6], 
in charm sector could be a Ac-A^ bound state with reasonable amount of binding energy, 
which deserves a further investigation. Our result also tells that the newly observed 
"exotic" state in bottom sector, the 10890) [37j, could be treated as the Ab-Ab bound 
state, whereas with an extremely large binding energy. 

It is worth emphasizing at this point that although our calculation result favors the 
existence of heavy baryonium, it is still hard to make a definite conclusion yet, especially 
with only the leading order two-pion-exchange potential. The potential sensitivity on cou- 
pling constants and energy cutoff also looks unusual and asks for further investigation. To 
be more closer to the truth, one needs to go beyond the leading order of accuracy in 
expansion; one should also investigate the potential while two baryon-like triquark clusters 
carry colors as proposed in the heavy baryonium model [HI [16]; last, but not least, the 
unknown and difficult to evaluate annihilation channel effect on the heavy baryonium po- 
tential should also be clarified, especially for heavy baryon-antibaryon interaction, which 
nevertheless could be phenomenologically parameterized so to reproduce known widths 
of some observed states. 
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Appendix 



In this Appendix, we present more detailed formulas and definitions used for the sake 
of reader's convenience. 

The 7 matrices take the following convention 

/ 1 \ / a' \ / 1 



And the Dirac spinors for read as 



where Xa is two-component Pauli spinor, and 



„(rt . , ™ f I . (64) 



2Ms V V. 



la 



where rja = —icr^Xa^ ^ind a = 1,2. Spin-| field for S"*"* is described by Rarita-Schwinger 
spinor u'^{p , a), which can be constructed by spin-1 vector and spin-i field [36], that is 



2My 



where is four-component Pauli spinor of a spin-| particle, and L(^)(p)(^ is the boost 
operator for spin-1 particle, 

( E Pj \ 

!r ^ (66) 

where i, j are indices of the space components of momentum p. The positive- and negative- 
energy projection operators for spin-| baryon are 

[^^{P)U = J2'''^{p,s)Mp^s) = (^^) (67) 

and 

[^~ip)]al3 = -Y^ Va{p, S)vp{p, s) = (^"^^T^^ , (68) 



±s 
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respectively. 

The positive- and negative-energy projection operators for spin-| baryon are 
[^Mp)] ccB ^a'- "(^^' *)^'^. ^(^^' 



±s 



and 

= - 5^ V/., a(p, S)V^, p{p, S) 



±s 



respectively. Here, and u are Lorentz indices; a and (3 are Dirac spinor indices. 
The basic Fourier transformation with Gaussian form factor reads 



Ar m\ „ f Kr m 



Svrr 

and hence 



e-'-erfc{-- + - -e-'erfcl- + ^ 



erfc{x) is complementary error function, which is defined as 

2 f°° 2 

erfc(x) — —p= / e~* . 
The factorization in double Fourier transformation goes hke 



H 
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(27r)^ ujiUJ2{uJi + a){uj2 + a)(a;i + UJ2) 



(27r)6 a^Vio K + A2)(a;i + A2) 

2 /""^ r/A 1 
Jo a + ^ a 
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Here, 



2a 

TT Jo (a 



)^ Ui{uJi + a) 



d^ki 2a 



^ki 2a r 

7r)3 TT Jo 



gikirg-k2/A2^j)^ 



(76) 



and for simplicity we define coi — V^i + '^2 = •\/k|~+m^ . 



In the case of the monopole form factor, i.e. /(k^) = 1 corresponding function 

to F(A, r) reads 



R{\ t) 



(27r)3 k' + m2 + A2 + k' + 



47rr 



(77) 



Operator Oi(ki, k2) contains two parts. The first part of Oi(ki, k2) while acting on 
functions in configuration space goes like 

Oi(ki, k2)F(A, ri)F(A, r2) = (ki • k2)2F(A, ri)F(A, r2) 

= (VhVi,)F(A, ri)(V2iV2,)F(A, r2) 

= ^F'(A, r)F'(A, r) + F"(A, r)F"(A, r) , (78) 

where 

V.V.^(.«-^)(i^)+^(|,). 

and the limit r2 — > ri = r is taken. The second part of 02(ki, k2) while acting on 
functions in configuration space goes like 

02(ki, k2)F(A, ri)F(A, r2) = (tri • ki x k2)(fr2 • ki x k2)F(A, n)F{\ r^) 

= (Tu(T2jeikiejmn{^lk^\m)F{\ ri){V 21V 2n)F{X, r2) 
^Ij^km^in ^Im^kn^ij ^In^im^kj) X 

(VuVi™)F(A, ri)(V2iV2n)i^(A, r2) 

_ 2 
3 



-4f'(A, r)F'(A, r) + -F'(A, r)F"(A, r) 



CTi ■ (72 



+ - 



2 /F'(A, r) 



- F"(A, r)^ 1f'(A, r)5i2 , 



(80) 
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where <Ti • a-i gives spin-spin potential and 5'i2 = ^^^2 — a\ - cr^ gives the tensor 
potential. 
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